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Abstract 

A new method of path averaging for waves propagating in a random dilute 
system of identical scatterers is developed. The scattering matrix of such a 
system is calculated. The method systematically takes into account repeating 
scatterings on the same scatterer and effects of correlations. Results obtained 
show the significant influence of new effects on the extinction of a coherent 
field and are valid in both the diffusive and localized regimes. 
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I. INTRODUCTION 



Diffusion and localization of waves propagating in random media have attracted a great 
deal of interest within the last decade [|]|§. Even though considerable progress in under- 
standing these phenomena has been achieved, many questions remain open. A commonly 
accepted picture is as follows. A coherent wave incident on a medium breaks into an inco- 
herent stream of photons within the distance of the order of the photon scattering mean-free 
path, i which then travel randomly through the medium. Thus, one can expect a smooth 
diffusive distribution of the wave amplitude and strong fluctuations of the phase at space- 
time scales of the order of the photon mean-free path and the free-flight time. When the 
wavelength, A, is much smaller then £, the transport equation or the random walk approach 
give an adequate description of the photon diffusion in a medium ||. On the other hand, 
when A and i are comparable, wave interference can destroy ordinary diffusion and can 
lead to the wave localization. To incorporate interference effects into the theory of wave 
propagation in random media one must start with the solution of the wave equation, 



Here ui is the wave frequency and the function £(x) describes properties of scatterers and 
their spatial distribution. When self- averaging of a field takes place in a random medium, 
observable quantities can be obtained upon averaging over an ensemble of random configu- 
rations of scatterers. An averaged Green function of Eq. (1), (G x ), defines the extinction 
length, £ ex , and the density of states. The transport properties of waves are described by the 



field-field correlation function, (G X G X ,). With the help of the Dyson equation, analysis of 

(G x ) is usually reduced to calculations of the mass-operator. The field- field correlator sat- 
isfies the Bethe-Salpiter equation. The latter involves the vertex function which describes an 
effective interaction of waves in a random medium. Microscopic analysis of these functions 
in a system of identical scatterers is based on multiple scattering perturbation theory. The 
perturbation procedure represents all quantities of interest as infinite series, where each term 
involves free-space Green functions, one-particle scattering operators and, partial distribu- 
tion functions of scatterers. Each term in the series corresponds to a particular sequence 
of scattering events and has a multiplier n m , where n is the average density of particles 
and m is the number of different scatterers participating in the scattering process. The the- 
ory contains three fundamental length scales, namely, the scatterer size, a, the separation 
between scatterers, n -1 / 3 , and the wavelength, A. In the case of dilute systems, when the 
volume fraction of scatterers, na 3 , is small, one can ignore all correlations in a scatterers 
distribution and use only the first order terms of the mass-operator and the vertex function 
as the leading terms of the expansions over small n. Most of the analytical results are based 
on this leading approximation in the low density limit and on simple corrections to it. 

The Green function and field-field correlator obtained in the leading approximation only 
takes a a restricted class of scattering diagrams into account. Only multiple scattering pro- 
cesses in which a wave scatters no more than once off any scatterer are included. This 
prescribes a "long memory" to a wave, such that when traveling in a medium it "remem- 
bers" all scatterers that were already visited and avoids them in its subsequent motion. 
In application to the wave diffusion, this restriction seems to contradict to the underlying 
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picture of randomly walking photons. In the case of particle-like photons this would lead to 
anomalous diffusion, significantly different from ordinary diffusion Besides, the fraction 
of these "long memory" paths is exponentially small for long path lengths. Therefore, the 
leading approximation must be accompanied by strong physical reasons for ignoring other 
statistically dominating scattering sequences. Such reasoning usually appeals to destructive 
interference of scattered waves which separates wave diffusion from diffusion of particles. 
This was analyzed in numerical calculations on finite clusters of scatterers M, but reliable 
analytical results in this area are absent. The problem arises because, in order to investigate 
an influence of wave interference, one must evaluate the entire series for the Green function 
and the field-field correlator. 

In this paper we propose a new method that enables us to consider processes with 
repeated scattering in calculation of the scattering operator. Our method is based upon 
(i) averaging over random paths instead of averaging over ensemble configurations, (ii) the 
statistical independence of random jumps of scattered waves, (Hi) the existence of a hierarchy 
of length scales in dilute systems of scatterers. 



II. LOW DENSITY LEADING APPROXIMATION 

The Green function of Eq. (1) corresponds to a field induced by a 5-source. It satisfies 
an integral equation with the iteration solution given by the formal operator series: 

/\ A A A A A A A A A A A A A 

G =D +D£D +D£D£D+ . . . =D +DTD , (2) 

where D is a free-space propagator and T denotes the total scattering operator of a system. 
One can choose any representation to obtain a constructive form of these expressions. For 

instance, in the coordinate representation, operators D and £ have the following matrix 
elements: 

e ig|x-x'| 

D XX ' = — : 77 , £ XX ' = g 2 £(x)5 xx / , (3) 

47r|x— x'| 

where q = u is a wave number in the system of units where speed of waves c = 1. 

In the case of a single scatterer, using an exact solution of Eq. (1), one can obtain the 

scattering operator of a single particle, t. In a system of TV identical scatterers, the total 

A 

scattering operator can be expressed in terms of t-operators 0: 

A OO A A A A A A 

T = E E t x D...D tDtx (4) 

s=l (x r ..x s ) 

A 

where t x denotes the scattering operator of a particle located at the point x m . The 
sums in Eq. (4) account for all scattering paths, (xj— > x 2 x s ), that are possible 

within a given spatial distribution of particles. Ensemble averaging of Eq. (4) eliminates 

this restriction and gives an average scattering operator (T) ensemble = T in a form of the 
weighted sum of contributions from all scattering paths: 
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T = E / ^(x 1 ...x s )t x D...Dt X2 Dt Xi , (5) 

where symbol J^ x x ) denotes integrations over the all indicated variables with the natural 
measure for each variable. The path weight, W(~x l . . . x a ), is the fraction of those configu- 
rations where the path is allowed and it can be expressed in terms of partial distribution 
functions of scatterers. 

The conventional way of calculating T in a statistically homogeneous system makes 
use of the momentum representation where the free-space propagator is diagonal, D kk , = 
5 kk , D£ = S kk , [k 2 — to 2 } , and the scattering matrix of a single particle depends upon the 

A 

location of a scatterer in a simple way (k'|t x |k) = t k , k exp ix(k'— k). Ignoring all short- 
range correlations arising from the rigid nature of scatterers, one can calculate exactly the 
"long memory" subseries. This subseries contains those terms of Eq. (5) where the scattered 
wave visits different scatterers only, so, all x m are independent and W^x-^.-xJ = n s . 
Retaining only these terms is known as the low density leading approximation or as the 
approximation of independent scattering. It is thought to be appropriate in the case of dilute 
systems of small scatterers and most of analytical results are based on it. For instance, the 
diagonal part of the T-operator in this approximation has the form H: 



rplcading _ 4- 



l-nt M D£ 



(6) 



where t>k,u — (k|t |k) is the scattering amplitude of a single particle in the forward direction. 

However, even a naive statistical analysis gives rise to doubts regarding this approxima- 
tion. Let us consider a system with the linear size L containing N scatterers. A wave 
propagates through the system by means of random jumps from one scatterer to another. 
During this motion a likely jump length is of the order of the mean-free path, i. In the 
diffusion limit after s jumps a wave covers the distance £^/s and leaves the system when 
s > Sq = L 2 £~ 2 . Therefore, paths with length of the order of So are most responsible for 
the diffusive content of an emerging field. A total number of paths consisting of s jumps is 
equal to N s . The number of "long memory" s-paths is N\/(N — s)\. It is straightforward 
to estimate that in a macroscopic system the ratio of these two numbers, or fraction of the 
"long memory" paths, tends to zero as exp — ^/s/2n\/J . Thus, all "long memory" paths 
represent a statistically negligible fraction of the scattering paths that are important for the 
wave diffusion. 

Moreover, the "long memory" in random walks provides an effect of "outward pressure" 
acting on a walker. For a simple walker the probability to stay within any volume around 
the origin after it makes s steps is proportional to N s /V s = n s . When a walker is 
forced to avoid repeating scattering on the same scatterer, this probability is proportional 
to N\/ (N — s)\s\V s ~ n s /s\. The longer paths we consider the smaller, comparing with 
an unrestricted walker, this probability is. In the latter case an average distance from 
the origin, R s , is of the order of i\fs. Since for a particle with the "long memory" the 
probability to stay within any volume around the origin is weakened by the factor 1/s!, 
the R s must be greater than that for a simple walker. This can be treated as an effective 
outward pressure leading a particle much farther away from the origin. The similar pressure, 
caused by excluded volume, exists in the case of the self-avoiding random walks and it leads 
to anomalous diffusion fH|. 
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From this point of view the success of the leading approximation in a description of the 
wave diffusion looks rather surprising. In any case, it is necessary to investigate the effect 
of all remaining paths on the scattering operator and the field-field correlator. 

III. RANDOM PATH AVERAGING 

With the use of uncorrelated distributions of scatterers one can calculate the "long 
memory" subseries exactly. However, it does not help in the evaluation of the total scattering 
operator without bias against any scattering paths. To advance here we modify the averaging 
procedure. Averaging over random paths of scattered waves is physically equivalent to 
averaging over random distributions of scatterers. Therefore, we can rewrite Eq. (5) using 
jump vectors, R m = x m+1 — x m , as independent variables instead of coordinates of scatterers 

T = E / R R W ... +R D. . . Dt x+R+R Dt x+R Dt XQ . (7) 

Any possible scattering path, {x , R x . . .R s }, is accounted in Eq. (7) with a path weight, 
W(xq, Rj. . .R s ). In the system of N scatterers fixed in space, the total number of s- 
paths is equal to (N — l) s and, as one can see from Eq. (4), they are all accounted for 

in T. When scatterers are randomly placed in space, a path weight is proportional to the 
probability to find all scatterers in proper places along a path. This probability is given by a 
joint distribution function of s scatterers which takes into account all correlations between 
them. If one ignores these correlations, a path weight takes a form used in the leading 
approximation, W = (^r") =n s . Any attempt to preserve all correlations in W presents 
too general problem with no constructive way to go far beyond the leading approximation. 
Meanwhile, a random path picture allows us to approach the problem utilizing the random 
jump distribution instead of distribution functions of scatterers. On each random step along 
the path {x ,R 1 . . .R s } a wave may jump to any of N—l scatterers. If we introduce a 
density of scatterers available for a single jump, ng(R), and assume statistical independence 
of random steps, then the path weight can be factorized as follows: 

W(x , R r RJ = n^giRj. . .g(R s ) . (8) 

The function ng(r) discribes the density of scatterers at a distance r from the one placed 
at the origin. In a statistically homogeneous system this function is close to the average 
density, n. However, correlations between scatterers and fluctuations of their concentration 
cause deviations of g(r) from unity at small distances. This is supported by the following 
argument. The fluctuations of density are discribed by the Poisson distribution. In a dilute 
system an average volume per particle, n _1 , is much greater than the volume of a scatterer, 
a 3 . Therefore, analyzing density fluctuations at large scales we can neglect the volume of 
the scatterers themself. The probability to find m particles within a spherical shell of 
radius r and thickness n" 1 / 3 is P(m, r) = (4irr 2 n3) m exp (— ^nr 2 n^) jm\. It shows that the 
probability of large fluctuations of density rapidly decreases at r> n -1 / 3 . We can expect, 
therefore, that the deviations of g(r) from unity are substantial at distanses comparable to 
the average separation between particles, n -1 / 3 . In application to wave propagation further 
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simplification seems to be relevant. When the scattering crossection, <7 SC , is small the precise 
structure of ng(r) more likely can be disregarded since a wave with a large probability 
can just ignore a lot of surrounding scatterers. The probability that a wave makes an r- 
jump is proportional to the probability of not meetting any scatterer inside a cylinder of 
lenght r with the base area equal to <x sc . From the Poisson distribution it follows that 
P(0, ra sc ) = exp (— nra sc ). Therefore, we can assume that the effect of density fluctuations 
on ng(r) is substantial at distances of the order of the scattering mean-free path, (na sc ) _1 . 
The function that takes into account the above arguments can be choosen in the form 



ng(r) = n 



7 3 

1 - - — exp(-r7) 
8nn 



(9) 



This function is normalized on Af-l and, consequently, the path weight W(xo,R . . . RJ 
is normalized on the total number of s-paths. A parameter 7 is understood to be of the 
order of na sc , if the scattering mean-free path is much greater than n -1 / 3 , and it tends to 
the inverse average separation, n^ 1 / 3 , when the scattering cross section becomes large and 
both lengths become comparable. 

In a dilute system Eq. (7) can be further simplified by noting that the size of a scatterer 
is much smaller than the average separation. If we use the coordinate representation in all 

intermediate states for the matrix elements (k'|T|k) = T k , k , 

T, k = E/ x Rl WM^.+ .+K,)-^ t k , r ,D r ,_ r .«....D r ,_ ri+R t rik , 

s=0 > / (x ,±ti...±t s j J (r 1 ,r 1 ...r s ,r s j 

(10) 

we can see that the ranges of changes of variables (Ri . . . R s ) and (r 1; r[ . . . r s , r'J have 
different order of magnitude. A single particle scattering matrix in the coordinate represen- 
tation, t r ' r , vanishes outside the scatterer, where r' r > a, while the distances between two 
successive scatterers, R, are of the order of the average separation, n" 1 / 3 . Assuming that 
short-range density fluctuations are not crucially important in the low density limit, we can 
assume that essential values of R are much greater than those of r', r. This fact allows 
one to use the far-zone asymptote of the free-space propagators in Eq. (10): 

iqR 

IWr' « V R e^'\ D R = _ , (11) 

where q = wn, n = H/R. It makes it possible to calculate all integrals with respect to 
( ri , r s , r'J in Eq. (10): 

T k , k -E / ^ k (x +R s +...+R 1 Hkx otk ^. . . t D . .-B Ri . (12) 

In a homogeneous system, the integral over the location of the first scatterer, x , gives 
the common multiplier <5 k / k expressing homogeneity of the system. Since the distribution 
function, g(R), depends upon the jump length only, the integrations over the directions of 
the jumps can be separated from the integrals over jump lengths in Eq. (12). All directional 
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integrals can be evaluated with the help of the Taylor expansion of smooth functions in the 
integrand and a deformation of the integration contour in the complex plane: 



(T=± 

This formula works quite well in the case of weak scattering anisotropy, 

at 



d(h • k) 



^4-<fci2t fi , (14) 



especially when /ci? ^> 1. Eq. (13) becomes the rigorous equality for isotropic scattering. 
After utilizing Eq. (13), all the remaining integrals over jump lengths in Eq. (12) can be 
factorized since there is no correlations beween jumps. The resulting expression for diagonal 
elements of T k / k takes the form 

oo 

T fcw =nt kk +n£ S t toT a 9v.K.a^ 1 ■ ■ ■ g^K^g^t^ ■ (15) 

8 = 1 °"l°"2'" 

In this equation the indexes a 1 ,a 2 ,--- correspond to "+" or "-" and we denote 

tkcr = ( k |t|o-a;k), t ffk = (<rwk|t|k) , 

-.977-77 poo 

t aa , = (<ju;k\t\<j'u;k) : g a = — j rg(r)V?e* kR °dr . (16) 

Due to the fact that all intermediate states in Eq. (15) are aligned along or against the 
incident momentum k, the multiple-scattering part of the T-matrix can be rewritten in 
terms of 2 x 2-matrices t and g defined in Eq. (16): 

00 00 

TL ultiple =n j- [g (tgrU , n £ f g (1-tg)- 1 ] ^ t CT , k . (17) 



s=0 00' s=0 a a' 



Finally, using an explicit expression for the matrix (1 — tg) , we can represent the T- 
matrix in the form 



l-tr||tg||+det||tg|| 



where tr ||tg|| and det ||tg|| denote the trace and the determinant, respectively. For a spher- 
ical scatterer t ++ = t is the forward scattering amplitude and t_ + = t + is the 

backscattering amplitude with all momenta confined at the "mass shell", k 2 = u; 2 , while in 
t k± and t ±k one of the momenta (k) can be of an arbitrary magnitude. 

This T-matrix is formally defined at all spatial scales and for any wavelength. In our 
evaluation, we assumed a weak anisotropy of scattering and a low density of particles, 
na 3 1. An isotropic scattering dominates, for instance, in the Rayleigh limit when a 
wavelength is much greater than the scatterer size, au 1. In this case, since scattering 
amplitudes do not depend upon directions of incoming and outgoing momenta, Eq. (18) can 
be simplified as follows: 
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' kui 



n 



4- nT l 



■i.9-\ 



0- 



(19) 



where t w = (o;k|t|o;k) is the on-shell scattering amplitude T ku = (k|t|u;k) is the on-off- 

shell amplitude and t fcw = (k|t|k) is the off-shell one. All these quantities can be obtained 
from the solution of the boundary-value problem for a single scatterer and are well-known 
for scatterers of a simple shape Quantities g± are given in Eq. (16) and depend upon 

the choice of the density of surrounding scatterers. The substitution of the function g(r) 
given by Eq. (9) yields the total scattering operator of a dilute system of identical isotropic 
scatterers: 



■ k u) ' 



n- 



^ku+ n 


T 2 ft 






U k Snn^k 





(20) 



This result involves the effective correlation radius 7 _1 . In the limit 7 — > 00, when 
correlations between scatterers can be ignored, T-matrix given by Eq. (20) reproduces the 
expression for the coherence length obtained in the leading approximation. The choice of 
7 ~ f" 1 physically means that the contribution to the coherent field from scattering paths 
consisting of many short jumps is assumed to be small, since they are rare in a statistically 
homogeneous dilute system. 



IV. EXTINCTION LENGTH 

The averaged Green function describes the coherent part of the field created by a point 
source placed in a random medium. Due to random scattering of propagating waves the 
pole of the Green function is shifted away from the pole of the free-space propagator, k = 
uj. The imaginary part of this pole defines the extinction length of coherent waves in the 
random medium. According to the relationship between Gfc, w and the scattering matrix, 
Gku) = D£ + (D£) 2 , it seems that G^^ has a pole of the second order when k tends 
to u. However, this singularity disappears, since T^^ has there a zero of the first order, 
T kljJ = —2u>(k—uj) + ((k — uj) 2 ). Our result, as well as ms , demonstrates this behavior. 
Therefore, only a pole in T-matrix gives an actual singularity of the Green function. In the 
leading approximation this pole can be found from the equation: 

k 2 - u 2 = nt k> „ (21) 

This approximation relies on the low density of particles and weak scattering, which implies 
that /= 4^a 3 n 1 and is small. Weak scattering takes place in the Rayleigh limit 

where aw C 1 and the scattering amplitude itself is a small quantity. An expansion over 
the small parameter a uo yields the leading term of the scattering amplitude 

Ret fcA) w £^7raV, (22) 

where e is the difference between dielectric constants of the scatterer and the host medium. 
The Ret fccJ does not depend upon k and is of the second order, while the leading term of 
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the imaginary part is of the fifth order, as it follows from the optical theorem, 



Im t fciW = (2tt 2 ) 2 cj J |t kk , 



dkf. 



Solving Eq. (21) for long waves in a dilute system we obtain 



k 2 {uj) = uj 



l + fe + i{fe)'- 



UJ 



Ann 



When /e C 1, this equation gives the well-known expression for the extinction length, 

= Im k{uj) = na sc /2^iuj(fe) 2 uj 3 /8nn. 

Here the factor uj 3 /8nn is of the order of (R/X) 3 , where R and A are the average separation 
between scatterers and the wave length, respectively. It can be dangerously large, but in the 
weak scattering regime it is totally neutralized by the small parameter fe. However, even 
in a rare mist of water drops this parameter can be of the order of unity due to the large 
value of e. In similar dilute systems with large mismatch indexes the leading approximation 
still seems to be valid. For R 3> A it gives 

l£~u>fe (R/\f^uj. 

On the other hand, the extinction length (the phase coherence length) being much smaller 
than the wave length is physically meaningless since then we cannot talk about wave prop- 
agation anymore. Therefore, the leading approximation requires that at least two condition 
to be met, namely, auj <C 1 and fe<^l. 

Our T-matrix was obtained under the conditions / < 1, aw C 1 and it involves a 
phenomenological parameter, 7 = £~^ TT . Our equation for the pole of the Green function 
reads: 



k —uj 



k 2 — (uj+i^Y 



■ n\ 



k 2 — (uj+i^Y 



r 



8nn 



(k 2 -u 2 ) 



(23) 



To be consistent, we have to use the Rayleigh scattering amplitude given by Eq. (22). When 
/£< 1, the right hand part of Eq. (23) is small and we obtain 



/-1 



7- 



For fe ~ 1 and R 3> A, taking into account that X/i COTT 1 and ^ 3 /8jcn ~ (R/i CO rr) 3 ~ 1, 
we have: 



I 



-1 



UJ" 



87m 



7 



inn 



In both cases we obtain physically meaningful expressions for the extinction length since we 
have assumed that 7 = ^^ approaches na sc in the weak scattering regime and it becomes 
proportional to n 1//3 when the scattering cross section becomes large. 
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V. CONCLUSIONS 



In summary, we developed a new method for studying interference of multiple scattered 
waves that takes into account multiple returns to the same scatterer. We evaluated the 
scattering matrix of a dilute system of identical particles. The T-matrix obtained in Eq. (20) 
is valid in the Rayleigh limit, uoa <C 1. We did not impose any other restrictions on the 
wave length, therefore, our result can be used to analyze the localized regime, iot ~ 1. All 
scattering paths and finite range of correlations are accounted for in our approach. This 
results in the momentum cut-off that appears in Eq. (20) as a truncation of the free-space 
propagator, D£, by the term g^-D£ +l7 . It also gives an extra term in the numerator of 
Eq. (20), which is missing in x lcading . In the Rayleigh limit this term is small because the 
leading terms of all scattering amplitudes are equal. However, the change of the denominator 
is physically important for the consistent description of the wave extinction in a system with 
large mismatch indexes. Using a microscopic expression for the density function, g(r), in Eq. 
(19) one can obtain an exact expression for the extinction length. According to the definition 
of g(r), this function must have the structure g(r) = 1 — "decay function". It guarantees 
the appearance of the momentum cut-off in the denominator of the exact T-matrix. This 
effect is due to scattering paths with multiple returns. Even if each such path gives a 
negligible contribution, they altogether present a dominant amount of scattering paths and 
they cannot be neglected in the coherent field. Referring to the motivation of our work 
given in the introduction, we can conclude that the interference cancellation of non-"long 
memory" contributions in the coherent field is overestimated in the leading approximation. 

The T-matrix itself does not allow one to analyze the influence of "long memory" effects 
on transport properties of waves. Calculations of the field-field correlator are currently in 
progress. Our method enables us to proceed without utilizing the Bethe-Salpeter equation 
and does not require any additional approximations. 

We wish to thank E. Kogan for the useful discussion and A. Z. Genack for reading and 
commenting on the manuscript. This work was supported by the NSF under grant No. 
DMR-9632789, by the CUNY collobarative grant and by PSC-CUNY research awards. 
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